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ABSTRAGT 

q'' We develop the canonical theory of gravitational collapse in 2+1 dimensions with a negative 
CT. cosmological constant and obtain exact solutions of the Wheeler-DeWitt equation regular- 
ized on a lattice. We employ these solutions to derive the Hawking radiation from black 
holes formed in all models of dust collapse. We obtain an (approximate) Planck spectrum 
near the horizon characterized by the Hawking temperature Th = h\/GKM /2tt, where M 
^ ', is the mass of a black hole that is presumed to form at the center of the collapsing matter 
cloud and —A is the cosmological constant. Our solutions to the Wheeler-DeWitt equation 



are exact, so we are able to reliably compute the greybody factors that result from going 
beyond the near horizon region. 
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•rH , I. INTRODUCTION 

X 

. 5t , A principal goal of modern theoretical physics is the construction of a consistent quantum 

theory of gravity. In fact much effort has been directed at this problem over several decades 
with the result that there are at present several proposals on the table. However, it is 
probably fair to say that none of the proposals currently under development is either fully 
understood, free of ambiguities or universally accepted [ll]. 

In the absence of either a full quantum theory of gravity or of any direct experimental in- 
put, it seems worthwhile to address the quantization of particular models by the application 
of as wide a variety of techniques as possible. Among the most interesting and best under- 
stood models are those with spherical symmetry, while an especially interesting model with 
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spherical symmetry is that of self-gravitating matter. One model of gravitational collapse 
that is very well understood on the classical level is that of Lemaitre 0], Tolman j^, and 
Bondi [3| (LTB), which describes a self-gravitating cloud of pressureless dust. The system 
is known to develop both covered and naked singularities and therefore a successful quanti- 
zation of these models has the potential to address also the fate of the naked singularity. 

Ideally, one would like to have a quantum model that is able to predict what the final state 
of collapse will be, in particular whether or not and under what conditions the formation of 
a singularity may be avoided ^ . The model should also explain some features of the end- 
states. For example, it has been known for a very long time from semi-classical arguments 
that black holes behave as heat reservoirs and evaporate as black bodies with a characteristic 
temperature that depends only upon their conserved charges 

Baa. 

Explaining these 

thermodynamic properties in terms of their microstates would be a desirable feature of 
any quantum theory of collapse. Further, semi-classical arguments indicate that the nature 
of particle creation from a naked singularity is significantly different from a black hole 
and possibly even that they are unstable j^; however, because the semi-classical treatment 
breaks down when curvatures reach the Planck scale, a full quantum treatment is essential 
to understanding their evolution. This is another area upon which a good quantum model 
of collapse should shed some light. 

In order to examine these problems we have set up a canonical quantization program 



Id, lUl to describe a self-gravitating dust cloud using the variables introduced in pL2]- In 
our work, the classical geometrodynamic constraints of the gravity-dust system were given in 
terms of a canonical chart consisting of the mass contained within spherical shells, the area 
radius, the dust proper time, and their conjugate momenta. The diffeomorphism constraint 
was used to eliminate the momentum conjugate to the mass function and this procedure was 
shown to result in a much simpler constraint that was able to take the place of the original 
Hamiltonian constraint. Dirac's quantization then led to the Wheeler-DeWitt equation for 
the wave functional describing the quantum collapse. Later it was shown that the WKB 
treatment of the Schwarzschild black hole in this canonical picture describes Hawking ra- 
diation [l3|. To go beyond the WKB approximation, the Wheeler-DeWitt equation was 
then regularized on a lattice and quantum corrections to the radiation were proposed in 



14| . However, the lattice regularization turned out not to be differomorphism invariant. 



This issue was addressed in [15[. When care is taken to ensure that the momentum con- 
straint is fulfilled in the continuum limit, the lattice wave functional becomes described 
by the Hamilton-Jacobi equation and two additional constraints which uniquely determine 
the factor ordering and make it possible to obtain exact solutions for the wave functional. 
Subsequently, we reexamined the Hawking radiation for the general case of non-marginal 



collapse [I6j, recovering the Planck spectrum near the horizon and additional grey-body 
factors as the near-horizon approximation was relaxed. 

In this paper we turn our attention to the situation in which a negative cosmological 
constant is present. We feel it is desirable to first adapt our program to the lower dimensional 
2-M-dimensional case, where dust collapse has been examined in detail on the classical 



171 . Il8l | and the semi-classical level 19|] by two of us. Moreover, the BTZ black hole 20 



21 



which is the unique classical end state of the collapse, is reasonably well understood on 
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the quantum level from the point of view of the AdS/CFT correspondence 22, l23(]. Thus 



it may be possible in the future to compare the degrees of freedom employed by the two 
approaches. Most of the calculations are in close analogy to the calculations presented in 
our earlier papers on the LTB model, but there are various subtle differences (boundary 
conditions, relative factors in expressions) that need to be addressed; therefore, we present 
the analysis in sufficient detail to provide a self-contained presentation. 

Our paper is organized as follows. In Sec. II we review the classical model of self- 
gravitating dust in 2+1 dimensions and set up the canonical formalism in ADM variables. 
In Sec. Ill we transform to the more transparent Kuchaf-type variables appropriate to 
the presence of a negative cosmological constant. Here we determine the mass function 
in terms of the ADM variables and go on to express the canonical constraints in a chart 
consisting of the radius, i?, the mass function, F, the dust proper time, r and their conjugate 
momenta. We also take special care to address the boundary action. In Sec. IV we obtain 
the Wheeler-DeWitt equation and present the exact quantum states appropriate to a lattice 
regularization. In Sec. V we show how Hawking radiation, together with the appropriate 
grey-body factors, can be recovered from our exact solutions. We summarize our results in 
Sec. VI and conclude with a few comments on potential future developments. 



II. THE MODEL 

A. The classical solutions 

We are concerned with a self-gravitating, pressureless dust cloud, described by the energy- 
momentum tensor T^;^ = e U^Uv, in 2+1 dimensions with a negative cosmological constant 
—A, A > 0. The metric may be given in comoving, synchronous coordinates as 

ds^ = -dr^ + e^"^^''^dp^ + Rir, pfdip^ , (1) 

where r is the dust proper time, and p labels dust shells of physical (curvature) radius 
R{t,p). Inserting this line element into Einstein's equations leads to the relation [l8| 

2b{T,p) _ {dpRf 

^ ~2{E-Fy ^> 

where E{p) and F{p) are time independent but otherwise arbitrary functions of p, which 
satisfy 

27rGe(r,p) = -^|^, {d^Rf = 2E - AR\ (3) 

We note that the gravitational constant G has the physical dimension of an inverse mass in 
2+1 dimensions (we set c = 1 throughout); the quantities F and E are dimensionless. 

The case of collapse is described by the condition drR < 0. There is still a freedom to 
rescale the shell label p: this can be fixed by demanding that 

R{0,p) = p, (4) 
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so that at the initial epoch (r = 0) the label coordinate is equal to the curvature radius 
R. This allows us to express the functions E{p) and F{p) in terms of the energy density at 
r = 0. From ([3]) we find 

F{p) = 2nG r p'e{Q,p')dp' 
Jo 

E{p) = [drR{0,p)f + Ap' (5) 

The physical interpretation of these relations is that 2F{p) represents the gravitating mass 
inside the shell labeled by p, and E{p)/2 is the total energy per unit mass of that shell. 
Therefore F{p) is generally known as the "mass function" of the shell and E{p) as its 
"energy function" . 
The solution to ([3]), 



R{r, P) = Y X '^"^ (^-^^ + sin-' y^py (6) 
shows that the shell labeled by p reaches a curvature radius i? = at the time 



At r = t{p) the shell becomes singular; therefore, r can take values no larger than tq{p). 
A detailed analysis, performed in [l8|, shows that only those shells that obey the relation 
2F > AR^ become trapped. Not only does this mean that, for a black hole to form, F > 
but also that a shell becomes trapped when it collapses to a size less than ^2F/A. 



The solutions can be matched to an exterior BTZ black hole 17, 18 



ds' = -{AR' - GM)dT' + + RV, (8) 

at some boundary that is specified by a fixed shell (the outermost shell), labeled pb- The 
matching requires the mass function at the boundary to be related to the mass parameter 
of the BTZ black hole according to F{pb) = GM/2. 



B. The Hamiltonian 

The ADM metric with circular symmetry for a 2+1 dimensional system takes the form 

ds'^ = -N{t,r)dt^ + L^{t,r){dt + N'drf + R^{t,r)d^'^, (9) 

where N is the lapse function and A^'' is the only shift that survives the symmetry. In terms 
of these metric components, the Einstein-Hilbert action can be written as 



= G-^ I dt I dr 



oo 



^{L'R' - LR" + AL^'R) - ^[{N'L)' - L][Ar^i?' - R] 

1j i V 



(10) 
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after integrating out the angular dependence. From here, we obtain the momenta conjugate 
to R(t,r) and L(t,r), 

Pr = ~lL-{N"Ly]. (11) 

and the action can be put in the form 

//*oo 
dt j dr[PLL + PnR - NW - N'H,] + (12) 

where Sqy, is a boundary term to be discussed below. Because no time derivatives of the 
lapse or shift occur in the action, they are Lagrange multipliers. The Hamiltonian and 
momentum (diffeomorphism) constraints are given respectively by 

---1 A D r I / ^ \ 



= -GPlPr - G-'ARL + j ^ 

Hr = R'Pr - LP'^ ~ 0. (13) 

The total action is the sum of (fT2l) and an action Sd that describes the dust. The canonical 
formalism for dust was developed in j2j| and elaborated in 25]. We consider only non- 
rotating dust, for which 

//•oo 
dt j dr[Prf - NW^ - N^W^], (14) 

where the Hamiltonian and momentum constraints are given by 



nt = r'Pr (15) 

The matter configuration will also act as a time-keeper for the quantum theory. In principle 
one could think to use a fundamental field for this purpose, but this would make the problem 
much less tractable and, as we will see, the main features of the theory are already contained 
in the dust model. 

It is easy to verify that the Poisson-bracket algebra of the constraints closes and the 
system is first class. With the total Hamiltonian denoted by H, the Hamiltonian equations 
of motion read 

R = {R, H}pB = 1^ = -NGPl + iV^/?', 
oPr 

L = {L,H}pB = ^ = -NGPr + (iV^L)', 
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5M / r'2 

i = {r,e}pB = — = N^l + — + iVV, 

Pn = {Ph, e}pB = = + ^ + G'-^A^A^ + iN'-PR)', 



AMP.H}PB^-f ^(^^p= + iV'P,). (16) 

and from them the equations in ([3]) can be recovered in the gauge t = t. 



The fall-off conditions 



We consider only mass functions which are such that at infinity the solutions approach 
the BTZ spacetime. This would be true in models in which the collapsing metric either 
asymptotically approaches or is smoothly matched to an exterior BTZ black hole at some 
boundary ph. It is then necessary to adopt the following fall-off conditions at spatial infinity: 



R{t,r)- 
L{t,r)- 

PL{t,r) 

N{t,r) - 
N''{t,r) 



GM+{t)r-^ 

AV2 ^ A372 

00/ -2\ 



A r - 

-2 



GM+(t)r-i 



A 



N+{t) + 0^{r- 



(17) 



On the other hand, the fall-off conditions at the origin depend sensitively on the conditions 
we place on the energy and mass functions. If we wish to avoid shell crossing singularities 
then the energy function must be positive, with negative slope. Near the center (r = 0) we 
take ^ = E„ ^n^"' with > and < 0. Likewise, take F = Y,^ F„r" with Fq > 0. In 
this way, we determine the following conditions as r — 0: 

R(t,r) a(t) + 6(t)r + C°(r3), 
L(t,r)^76(t) + O0(r), 
Pj,(t,r)^Pnoit) + 0\r), 
PL{t,r)-.PLo{t) + 0\T), 
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N{t,r)^^No{t) + 0^{r'), 

iV'^(t,r) ^ C°(r). (18) 

They are consistent with the constraints (fT3l) and are preserved by the equations of motion. 
With these fall-off conditions, the only non-vanishing boundary variations that arise read 

dtN+{t)M+{t), (19) 

where iV+(t) is the lapse function as r — oo and GM^ = 2F{r oo) is the ADM mass and 
a similar contribution at r = 0, 



j dtNo{t)Mo{t) (20) 



where A'o(t) is the lapse function as r ^ and GMq = 2F{0). To avoid the conclusion that 
and Nq freeze the evolution at the respective boundaries, the boundary terms must be 
cancelled by an appropriate boundary action. This can be achieved by adding the surface 
action 

^aE = - y dtN4t)M+{t) + j dtNo{t)Mo{t) (21) 

In the next section we show how this surface action can be absorbed into the hypersurface 
action by introducing a new canonical chart. 



III. CANONICAL TRANSFORMATIONS 

A. Mass function in terms of the canonical variables 

Now consider an embedding of the ADM metric in the spacetime described by the 
metric ([T]), 

ds' = -dr' + ^^f_p^ dp' + R'd^\ (22) 

where R = R{t,p), R = dpR, and we will use R* = drR as opposed to a prime and a dot 
for derivatives with respect to the ADM labels r and t, respectively. We shall set G = 1 in 

.et ] 

' 2 ^2 ,2 12 



the following to simplify the expressions. Let R = R/ a/2(-E — F), then 

= R%'^ - r" 



N=^{fp'-pT'), 

= ^P^^. (23) 



Inserting these into the momenta (fTTl) . we find that 

^ IkiR'p'^ -t'^)-(r'pp' -tt')R! . (24) 



LPl = -= 

R{fp' - pr') 



Derivatives with respect to the ADM time can be exchanged for derivatives with respect to 
the proper time using 



R 



R*f + Rp = R*f + Ry/2{E - F) p, 



R' = R*t' + Rp' = R*t' + R^2{E - F) p', 



which gives 



LPr 



R'R* 



(i?*2_2(E-F) , 
H , . = — T . 



V^iE-F) ^2{E - F) 
Substituting Einstein's equation ([3]), R*^ = 2E — AR^, then yields the simphfied form 

RW2E-AR'^ Ai?2-2F , 



LPl = T- 



V^iE-F) y/2iE - F) 



(25) 
(26) 

(27) 



which may be solved for r': 



LPlJ2(E-F) R'V2E - Ai?2 
T = . 1„ ^„ — - =F 



(28) 



Ai?2_2F ' Ai?2_2F ■ 
Substituting this into the expression for L in fl23l) . 

L2 = R%'^ - r'\ (29) 

and solving for F gives an expression for the mass function in terms of the canonical variables, 



F =- 
2 



Pl-^+AR' 



For future reference we introduce the function T defined by 
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T = AR^-2F = ?-- Pi 



(30) 



(31) 



JF = determines the apparent horizon and JF will play an important role in the quantum 
theory. One can check that though appears in the denominator of (125]) . r' continues to 
be well behaved across the horizon, as expected. (We note the difference in the definition of 
JF compared to the LTB case where one has JF = 1 — F/R (l5|.) 



B. New variables 



As in the case of the Schwarzschild black hole 12] and the LTB collapse model in 3+1 
dimensions lO], one can make a canonical transformation that elevates the mass function 
to a canonical variable. Interest ingly, t he expressions are similar to those presented by us in 
earlier papers on LTB collapse 15|, ll6|] . By directly taking Poisson brackets the momentum 
conjugate to the mass function is found to be 

LPl LPl 



AR^ -2F T 



(32) 
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Looking for a canonical transformation that would take the diffeomorphism constraint to 

Hr = R'Pr - LP'^ R'Pr + F'Pf, (33) 
we find a simple expression for P^: 

(34) 

where A = {LPl)'R' — {LPl)R" . Our momentum constraint becomes then indeed 

Hr = R'Pr + F'Pf ^ 0, (35) 

but we must first show that the transformation from the old set {R, L, Pr, P^} to the set 
{R, F, Pr, Pp} is canonical. We can do this by explicitly constructing the generator of the 
transformation. Denote it by F. Since we already know two coordinates {R and F) and one 
conjugate momentum (Pp), we use 

dF{r') 6¥ 



— f dF(r') S¥ 



0= Ur'PAr')WT\+ * 



0. / 



dPAr) 5PAr) 



The last equation in (136|) tells us that F = F[i?, L, Pl\. The third equation gives 



(36) 



^'^^'^ + 5PZ - ° ^ 5PZ - P|-P'VL2' 



and therefore 



whereas, from the first equation, 



¥= dr 



LPl - R' tanh- 



-1 LPl 
R' 



+ Fi[L,P], (38) 



showing that Fi = Fi[P]. Take Fi to be independent of P, a constant, and let us calculate 
Pr from the resulting F: 



Integrating by parts. 



,PfR'\' 5F 
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P , ARLPl (LPl/RJ 



T 1 - {LPl/R'Y 
_ ARLP, A 

we find precisely the new candidate momentum in ([MD- The transformation from 
{R, L, Pr, Pl} to {R, F, Pr, Pf} is generated by F. 

We now want to write the Hamiltonian in terms of the new variables, {R, Pr, F, Pp}- To 
do so, we use the two equations 



p/2 

P^- — - -T 



LPl = TPf. (42) 

Inserting the second into the first gives 

p2 -r-2 _ D/2 D/2 _ p2 -772 



and therefore 



as well as 



L2 ^ 



p2 ttS 

^ _ p2^2 ^^^^ 

Pr = Pr + ^RPf- p,^^_^,, . (45) 



Putting it all together, we find 



= -y [TPfPr + T-^R'F'] (46) 

which is of the same form as the expressions for the Schwarzschild black hole and LTB 
collapse in 3+1 dimensions, cf. Equation (42) in 15|]. The action in the new canonical 
variables then reads 

^-EH = j {PrT + PrR + PpF - NW - N'^Hr^ + Sq^ (47) 

with the new constraints (146|) and 

Hr = R'Pr + F'Pf. (48) 
We shall now discuss the boundary action in more detail. 



C. Boundary action 

The following considerations are in analogy to Sec. II. D in [151 . Because varying 
would lead to zero ADM mass and varying A^o would restrict Fq to zero, both and Nq 
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should be considered as prescribed functions. By the fall-off conditions, the lapse function, 
A^^, is required to vanish both at the center as well as at infinity. This implies that the time 
evolution is generated along the world lines of observers with r = constant. If we introduce 
the proper time of these observers as a new variable, we can express the lapse function in 
the form N^{t) = r+ and iVo(t) = Tq. This leads to 

^as = - j dtM+{t)T+ + j dtMoit)To. (49) 

Thus we remove the need to fix the lapse function at the boundaries. Extending the treat- 
ment in 121], the aim is to cast the homogeneous part of the action into Liouville form, 
finding a transformation to new variables that absorb the boundary terms. This can be 
done by introducing the mass density F = F' as a new canonical variable. Define 

F{r)=Mo+ dr'T{r), T{r) = F'{r), (50) 
Jo 

and reconsider the Liouville form 

POO 

e := / drPpSF - M+6t+ + MqSto 







POO 

= / Pf6F + T+5M+-To6Mo, (51) 
Jo 

where we have dropped an exact form. But 

SF = 5Mo+ [ dr'ST{r') (52) 
Jo 



gives 



(poo \ POO pr 

J dr'PF{r')-Toj5Mo + J drPpir) J dr'5T{r') +T+5M+. (53) 



Noting further that 
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POO pr POO POO 

/ drPpir) / dr'5T{r') = dr 5T{r) / dr'Ppir'), (54) 

Jo Jo Jo Jr 



^ See Kuchaf '12|. Consider 

/ 

dr'Sr{r') X / dr'PF{r') 



= 5T{r)x j dr'PF{r')~ Ppir) X j dr'Sr{r') 

Integrating the left hand side from to cxd gives zero, therefore 

drPpir) I dr'&V{r') = dr 6T{r) / dr'Ppir') 
Jo Jo 



11 



we can write the Liouville form as 



e = ^ dr'PF{r')-T^5M^ 

/oo / poo i-r \ 

dr6T{r) (j drPpir) - J dr'Ppir')] +T+5M+ 

poo noc nr 

+ {5M+-6Mo) / drPpir)- / dr5T{r) I dr Ppir') +T+1 
Jo Jo Jo 

POO 

= poSMo + + / Pr{r)5T{r), (55) 

where 

Po = -To, 

poo 

p+ = T++ / drPplr), 
Jo 

Pr(r) = - / dr'Ppir'). (56) 







The new form of the action is then 

^EH = j dt {poMo + P+M+ + j dr [PrT + PrR + PvV - NW - iV^H,]^ , (57) 
where the new constraints read 



= -\ [-^P^Pr + r'^B!T] + p^y'i + ^ ^ 0, 

Hr = R'Pr - TP^ + r'Pr ^ 0. (58) 

The Hamiltonian constraint can be simphfied if the momentum constraint is used to ehmi- 
nate Pr, and the constraints in (l58ll can be replaced by the following equivalent set 

P^ + :FPl-y^o 

Hr = r'Pr + R'Tr - ~ 0. (59) 

(Note that the Hamiltonian constraint has been 'squared' in order to arrive at the new form; 
it is therefore of dimension mass over length squared. Re-inserting G would correspond to 
the substitution F — T /G.) This involves a little algebra which has been described in the 
appendix of [l^. These equations will be used for quantization in the next section. 

Here we emphasize that the relative sign between the dust and gravitational kinetic terms 
can change because T is greater than zero outside the horizon and less than zero inside. This 
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change of sign is already present in (l58l) and has not been introduced by using the momentum 
constraint to ehminate Pp- It is of fundamental interest in the quantum theory because the 
Wheeler-DeWitt equation becomes locally elliptic outside and hyperbolic inside the horizon, 
which is of importance for the formulation of the proper boundary value problem. A change 



of sign was also noted for the case of a non-minimally coupled scalar field |2 
for an earlier discussion in the black-hole case. 



see also 



27 



D. Relation between dust proper time and Killing time 



In what follows, we will have need of the relationship between the dust proper time and 
Killing time; therefore we now address this issue. This subsection corresponds to subsection 



II. E in the LTB case 15|. We had obtained the expression for r' in terms of the canonical 
variables in Sec. IIA as 



R' 



F)Ty^2{E~F) 



T 



Defining a = Ij^liE - F) gives 



r 



a aj- 



(The equation of motion guarantees that the quantity \J\ — a^T is real.) 
value of a, this equation can be integrated. When R > a/2F/A, 

dR 



(60) 
(61) 

For a constant 



ar 



and 



ar 



A 



tan 



+ 



1 



aRy/A 



VI - V2AF 



: tanh" 



.1 yi - a^T 
R^K/2F 



A 



tan 



aRy/A 



+ 



1 



tanh 



_i R,/A/2F 



(62) 



(63) 



VI - o?T V2AF VT 

when R < a/2-F/A. We know that the spacetime surrounding a collapsing cloud in 2+1 
dimensions with negative cosmological constant is the BTZ black hole. The choice of sign 
depends on whether one is interested in an expanding (positive sign) or collapsing (negative 
sign) cloud. Matching the collapse and BTZ metrics also shows that T = a.t the boundary. 

Equations ( |62i) and (1631) may still be used so long as E' and F' are sufficiently small, since 
then we have a small amount of dust propagating in the BTZ background. In that case, it 
would give the relationship between the time used by families of freely falling observers and 
the Killing time, each family being characterized by a fixed value of E. 



E. Interpretation of the canonical data 



We have already reconstructed the mass function from the canonical data. In this section 
we want to reconstruct the energy function E and singularity curve Tq from the same. The 
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three functions E, F, and tq determine the collapse model completely, so reconstructing 
these quantities gives physical meaning to the canonical variables. 
If we begin with the momentum constraint in the form 

^'^^+^ (64) 



Pr P. 

and use the Hamiltonian constraint to eliminate P^, we get 



Substituting Pp = — Pp = LP^/JF, we find an expression, 



(65) 



that may be directly compared with the expression we had in (p8l) for r'. We see that 

Pr = , ^ (67) 
v/2(i?-F) 

and therefore the energy function is related to the canonical variables by 

E = =^ + F (68) 

2(rv^-^p|) 

where we have used the Hamiltonian constraint to give the result in terms of gravitational 
phase space variables only. Finally, knowing E and F in terms of the canonical variables, it 



is a simple matter to do the same for the singularity curve using the solution ([7]), cf. 15 



F. Hamilton equations of motion 

We now give the Hamilton equations of motion for the new system and derive Einstein's 
equations ([3]) from them. Introducing the smeared constraints 

poo 

H[N] = / drN{r)n{r), 
Jo 

POO 

Hr[N']= drN'{r),nrir) (69) 



and the Hamiltonian EI = T-[[N] + T-[r[N^], the canonical equations of the evolution are^ 
Mo = — = 0, 



6 Note that 

sr{r) ^ 

where 8 is the step function. 
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Po = 
M+ 

P+ -- 
f = 

Pr-- 
R = 

7r 



5W 

6m 



dr N ( pI + J, 



6p^ 



0, 



6M 



6M^ 



5m 



2NPr + N^'t' 



5Pr 

or 



5Pr 

6m 



— = 2iV kRP^ + 



sm 



5Pr 

p _ m 



(ATT)', 

/■r 

+2 / c/r'A^(r') 



To obtain Einstein's equations, we notice that because 



/\2 



N'P^. 



J'PrJ2{E - F) 



R N'^O 

f i 

the momentum Pr conjugate to R may be expressed as 

VR* 



(70) 



(71) 



P 



R 



(72) 

Note that we 



r^2{E-F)' 

Substituting this result in the Hamiltonian constraint leads directly to (E 
did not have to specify the lapse function. 

The lapse no longer has the interpretation of being the ratio between the proper time 
to BTZ time. The reason is that we have squared the original version of the Hamiltonian 
constraint. If we define a new version of the constraint by taking the square root, 



n, = Pr 



—2 r2 



then 



(73) 
(74) 



{T,n^[N^} = N^ 

shows that it is that carries that interpretation. 

Finally, we remark that the algebra of the constraints cannot be of the general form (given 
for example in [l|), again because we have used the momentum constraint to eliminate Pp 
in the Hamiltonian constraint. In fact, a short calculation gives 

{n[N],n[M]} = o, (75) 

{Hr [N^-] ,n[N]} = n [N^.N'- - , (76) 
{K[iVl, H,[Mn} = Hr [[N^, M^] . (77) 
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The Poisson bracket of the Hamiltonian with itself vanishes, in contrast to the general case 
where it closes on the momentum constraint. The other brackets coincide with the general 
case. The transformations generated by the Hamiltonian constraint can thus no longer 
be interpreted as hypersurface deformations. They are in general not orthogonal to the 
hypersurfaces, but act along the dust flow lines. 



IV. EXACT QUANTUM STATES 



A. Quantization and ansatz 



We now apply Dirac's quantization procedure to turn the classical constraints into quan- 
tum operators 

Px = -^T^, (78) 



'5X{ry 

and acting the operators on wave-functionals. The Hamiltonian constraint then reads 



vl/[r,i?,r]=0. 



(79) 



where A and B are smooth functions of R and F that encapsulate the factor ordering 
ambiguities. The factor ordering problem is unsolved and can be dealt with only after a 
suitable regularization procedure is implemented. We have included the formal expression 
5(0) to indicate the need for this regularization. Quantizing the momentum constraint using 



gives 



r 



5_ 



^[r,i?,r] = o, 



(80) 



but, as noted, the quantum constraints in (!79l) and (IHOj) are only formal until a regularization 
procedure has been selected. 



To regularize, we follow the construction 15(] and consider a one dimensional lattice of 
discrete points separated by a distance a which must be later taken to zero to achieve 
the continuum limit. We will require the wave-functional to (i) automatically satisfy the 
momentum constraint in the continuum limit and (ii) be factorizable into different functions 
for each lattice point. Such a wave functional is of the form 



^[r,i?,r] = exp 



I / drT{r)W{T{r),R{r),F{r)) 



lim TT exp [zcrri>Vi(ri, Ri, Fi) 



(81) 



where Fi = Yl]=o'^^i- Details of the construction are given in 15|]. One finds that the 
Hamiltonian constraint is satisfied independently of the choice of a only if the following 
three equations are simultaneously obeyed: 



aw. 



dRj J 



0, 
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± _|_ TT £ J- A ■ ■' 



' OR, 



0, 



5, = 0. 



(82) 

The first equation is the Hamilton- Jacobi equation and was used in [l^. The second equa- 
tion presents an additional restriction on solutions and the last equation tells us that working 
on the lattice is only possible if the factor ordering does not contribute to the potential term. 
If we find solutions to all three equations, we can do all other calculations on the lattice, since 
these solutions have a well defined continuum limit and satisfy the momentum constraint. 
We will now find exact solutions to all three equations. 

B. The measure 



The function A is closely tied to the inner product on the Hilbert space via the hermiticity 
of the operator 

It is therefore intimately connected to the choice of measure, fij. Below we consider the 

relationship between n and A, assuming that /i, is independent of r. For O = J^P^ to be 
hermitian, we require 



dRfx{R)^*{R){Oip{R)} 



dRfi{R){0^{R)}*iP{R), 



(84) 



'0 ^0 

which leads, provided the boundary conditions are trivial, to the (coupled) system 

2 d 



dR] 



d 



{^^,A,) = 0. 



(85) 



The second equation can be integrated to give 



1 d 



(/i,^,) - = A 



■J' 



and the first requires that 3?e(Q) = 0. We can view (15B]) as a (Bernoulli) equation for fij-. 



(86) 



1 _ Q,(F,: 

dR, \dR, V ' 



(87) 



whose general solution is 



±Q,(F,) / dR, 



+ a, 



(88) 
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where the upper sign is for the exterior (JF > 0), the lower sign is for the interior and 

l_(d\^ 



again with the same sign conventions. 

For solutions of (182|) . discussed in the following subsection, the relation between A and 

A, = \r,\dRMN\^3\)^ (90) 
which re-expresses ( l86l) with Qj{Fj) = (because Qj{Fj) can only be imaginary), will be of 
greater interest than (IHHI) . 



C. Exact solutions 

We have already emphasized that the signature in the kinetic part of the Hamiltonian 
constraint can change from elliptic (outside the horizon) to hyperbolic (inside the horizon) . 



This carries over to the kinetic term of the Wheeler-DeWitt equation. As discussed in [27 
we can say that the part inside the horizon is always classically allowed, whereas this is not 
necessarily the case for the outside part. The usual initial value problem appropriate for 
hyperbolic equations can thus only be applied for the region corresponding to the black hole 
interior. 

We determine the unique, exact solutions to (l82l) in Appendix A and here we summarize 
the results. Both in the exterior, > 0, and interior, JF < 0, we obtain a two parameter 
family of solutions: for the measure 

/i, = , (91) 

and for Wj 



Wj = const. ± ajTj ± / dRj^ , (92) 



where aj and jSj are constants of integration. The physical significance of aj can be demon- 
strated by acting on our wave functional with the dust energy operator P^: 



showing that a = 1/ y'2{E — F). The integral appearing in ( l92l) will be recognized then 
as the integral we evaluated earlier in (l63i) in order to obtain the relationship between the 
dust proper time and the BTZ time. The wave functionals are therefore oscillatory in the 
classically allowed regions. 

Given the requirement of factorizability on the lattice, there are no other solutions (either 
separating or non-separating) that solve the constraints, that is, we have obtained the com- 
plete class of factorizable solutions to all the constraints. Other solutions would necessarily 
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couple the shells composing the collapsing cloud. The solutions we have obtained, which are 
strictly valid on the lattice, contain two free parameters, aj and 6j, and we write them as 



y^. = e "''^^^e ' ' ^ ' \ (94) 

We know that a.j is connected with the energy and mass functions and that in general both E 
and F are functions of the radial label coordinate. It would therefore be natural to demand 
that a is also an arbitrary function of r, but this explicit dependence on r would violate the 
momentum constraint. The way out is to require the r dependence of a to appear via the 
mass function, -F(r), i.e., we must require that a = a{F{r)), and likewise b = b{F{r)). Any 
dependence on r via the mass function is allowed by the momentum constraint and since 
the Hamiltonian constraint does not contain a derivative with respect to F or F it continues 
to be obeyed. Hence, in the continuum limit we arrive at 



^[r, R, T] = e-/o°°^-r('-)b(i^M) exp J ±i drV{r) a{F{r))r ± j 



(95) 

The fact that a(r) is constrained by diffeomorphism invariance to depend on r only via F{r) 
means that the energy function is not arbitrary in the quantum theory but only those energy 
functions expressible in terms of F, i.e., E{r) = F{r) + l/2a^(F(r)) are allowed. 



V. HAWKING RADIATION 



A. Introduction 



The exact quantum states in ( !95l) describe the generic situation. In order to describe 
Hawking radiation, we need to introduce into the formalism the black hole and the analogue 
of the quantum fields which are used in the standard treatment. 

Following [l3| and [l6|, we consider a BTZ black hole surrounded by tenuous dust, by 
taking the mass function to be of the form 2F{r) = GMQ{r) + 2/(r), where 0(r) is the 
Heaviside function and 2/(r) represents the dust perturbation. Formally, the variable F = 
GM6{r) + f'{r) describes this perturbation and f'{r) describes the mass density outside 
r = 0. As described in 1^, the black hole state factors into a product state, one member of 
the which represents the black hole, while the other represents the dust and describes the 
Hawking radiation. 



The central idea of 13|] was to use the Bogoliubov transformation of the field operators in 
the BTZ spacetime. Here, since we have the exact quantum states at our disposal, we first 
identify the states that correspond to the ingoing and outgoing modes, respectively, of the 
standard approach and then calculate their inner product. The inner product is evaluated 
on hypersurfaces of constant Killing time, which correspond to observers who are static with 
respect to the black hole, not constant dust time, which would correspond to observers who 
are freely falling. We will therefore express our solutions in terms of the BTZ time, T, using 
([63]). 
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The relevant exact solutions of the Wheeler-DeWitt equation read (in this section we 
re-introduce h and G) 



± 



cir r 



ar ± 



(96) 



We will not concern ourselves with the normalizations at this point because they make no 
difference to our final result. 

For the positive frequency, incoming wave functional we choose from (jHS]) and replace 
the dust proper time with the Killing time using ( l63l) . We find 



(97) 



which has the standard form of a positive-frequency wave function. This wave-functional is 
independent of R. The lattice version of this state is obtained through the replacement of 
r = F' by a frequency uj, and is given by 



lim I I I 

o-^O J- J. 



X—o-bujj—icrujjTj)/Gh 



(98) 



For the out-going modes of negative frequency we have to take the state \E' . Inserting (!63|) 
into the corresponding state of (p6!l then gives 



T-2 I dR 



the corresponding lattice version being 



Tj-2 dR, 



1 - a]J^j 



(99) 



(100) 



B. The Bogoliubov coefficient in the near horizon Umit 



The Bogoliubov coefficients ly, |28| are given by 



13. 



Gh 



poo 



:ioii 



Using the coordinate transformation from (-R, t) to (/?, T), we have in {R,T) coordinates 

1 



RR 



a(Ai?2 - 2F) 



(102) 



where a 



, ^ Inserting the wave functionals (!98|) and fllOOp into (1101 1) then gives 

y/2{E-F) 



2uja 

'Gh 



^2{E - F) exp 



ah{uj + uj') ,aT{uj + uj') 



dR 



A/?2 _ 2F 



exp 



Gh 



2iauj , VT^a^' 
dR 



Gh 



X 



Gh 



(103) 
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We now change the variable of integration to T ^ keeping in mind the near horizon ap- 
proximation. T goes to zero on the horizon and diverges as R goes to infinity. We have 
dT jlKR = dR. Rewriting the integral in the new variable yields 



_2Lua 



dT = exp 

1T^k[T + 2F) 



2iaoj 




2j^^A(JP- + 2F)^ 

In the near horizon limit, we expand about !F — Q. The integral becomes 



(104) 



ah{uj + ijj') .aT{ijj + uj') 



dT T GhV2KF exp 



Gh Gh 

aoj fa?' 1 , ^ 



(105) 



GhV2AF V 2 4F 

Using the following formula after inserting the regularization factor e~^'^ for convergence, 

/ dx x'^-^e-(^+'«)^ = r(l/)(p2 + ^2)-i^/2g-i,.axctan(g/p)^ 

Jo 

^here ^ = ^ and g = ^ + ^) ^, we find 

GhV^KF ' Gh Gh J 

GhV2KFj \GhV2KF\2 AFJ) ^ ^ 

The absolute square of the above expression is given by 



and determines the particle creation via 

/•oo 

(in|iVout|in) = / d{aJ) |/3,,f . (108) 

(We integrate here over aio' in order to obtain a dimensionless expression.) Performing the 
integration and replacing ulo by GAe, where Ae denotes the energy of a shell, we get (setting 
also 2F = GM) 

, ,. , 7rAe(2E - GM)) e'^^^^/^ ^^^^^ 
(in A^out in) = ^ , " ■ ^= , (109 

where iVout represents the "out" particle number. 

We recognize that the Planck spectrum is modified by greybody factors which explicitly 
depend on the energy. These greybody factors are different from those obtained by taking 
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into account back-scattering (for example, see [29j) because we are using exact solutions of 
the quantum constraints. Both the gravitational field and the dust are quantized in our 
approach. 

The constant b that appears in the final expression for (in|iVout|in) also has its origin in 
the full quantum gravitational state. The simplest choice for 6 is 6 = 0, in which case the 
integral over uj' diverges. This divergence, however, is well-known and is connected with 
the normalization of the continuous modes. Dividing out the infinite constant from this 
integration the result, for 6 = 0, is^ 

m A/^aut m) = ^ , ; (110) 

This is a thermal spectrum with temperature given by 



/cbTh = n — - — , (111) 



which holds for each shell separately. 



C. Exact Bogoliubov coefficient 

It is possible to obtain an exact expression for the Bogoliubov coefficient as well by 
considering a series expansion about the horizon. We expand equation (11041) around JF = 
to all orders. Expanding the measure of equation (I104p . we get. 



Similarly the exponent can be rewritten as. 



n=-l 



def 

m=l 



GhV2AF 



ln^+5^/3^^"^ , (113) 



The integral can be performed before summing over the series of (I112p . 



'2{E -F) X 



GhV2AF 

n=-l 



^ We divide by Gh in order to keep the result dimensionless. 
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where /3i = -(l/a^ + 1/4F). 

m=2 



e (GHV2AF) — — = 1 + ^ JF"* 



-\ / ■ \ 2 oo oo 

^Hgot) E/".A^"^" - E7.^. (115) 



m,n=2 m=0 

where 70 = 1 and 71 = 0. The integral is 



00 «cx) 

Ghv2AF Jo 



:f/(g?i\/2AF) 

n,m=0 

V^iE - F) X 



GhV2KF 

n+m poo 



n,m=0 

and the exact Bogohubov coefficient can be written as 
GhV2AF 



(116) 

X=aujf}i/{GhV2AF) 



au ^— — / aT(uj + uj')\ 



(n \ -i(Tu}/{GHV2AF) / 
GhV2KFj \GhV2AFj 

(•_^)n+m I la^^ I X ... X ( '^^ -n-m+l] x 



X [1 + 



n+m=l 



f acoP^ Y 
\GhV2AFj 



(117) 



The simplest correction term to the near-horizon approximation is (^-independent and ob- 
tained for n = 1, m = 0: 

ao — — 1/4F and /3i = — (a^/2 + 1/4F), which can be simplified to be —j^- 

In general the higher order corrections are cu dependent and therefore lead to non-trivial 
modifications of the greybody factors. However, in the limit as a; — > 00, the corrections once 

again become independent of uj. These correction terms cannot be obtained in the standard 
derivation of the Hawking radiation because the geometric optics approximation is strictly 
assumed. We can obtain them because we have the exact quantum state at our disposal. 



VI. DISCUSSION 

In this paper we have extended the canonical formalism originally developed for the 
Schwarzschild black hole and later for the LTB collapse of inhomogeneous dust in 3-1-1 
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dimensions to 2+1 dimensional collapse with a negative cosmological constant. In this work, 
our emphasis has been on the physical parameters E, F, and tq, and we have succeeded 
in expressing them in terms of the canonical variables. This renders the formalism more 
transparent. 

Using a lattice regularization scheme that correctly implements the diffeomorphism con- 
straint in the continuum limit and assuming factorizability on the lattice, we presented exact 
and unique solutions to the quantum constraints and extended them to the continuum. They 
coincide with the WKB approximation, which therefore become exact. 

Using the canonical formulation we produced a general relationship between the dust 
proper time and the Killing time and then used our solutions to derive the Hawking radiation 
from a BTZ black hole. It was possible to reliably calculate the greybody factors because our 
solutions are exact. No further corrections to the Hawking radiation are obtained because 
the WKB solution is exact within the framework of our regularization. 

It is worthwhile comparing our derivation of Hawking radiation with the corresponding 



calculation in the four- dimensional LTB case 161, and also with the semiclassical 2+1-d 



derivation of Hawking radiation 19|]. There are significant differences in the gravitational 
collapse between the 2+1-dimensional and the 3+1-dimensional case on the classical level. 
The effects are indeed captured in the quantum formalism developed in the paper. Through 
canonical transformations both Wheeler-DeWitt equations can be transformed to the same 
form. The qualitative difference between the two cases is brought out by the difference in the 
form of JF which captures the information about the details of trapped surface formation. In 
the calculation of Hawking radiation, the framework is the same as in the four- dimensional 
case. In both cases it yields the correct Hawking temperature, though the correction terms 
for the Bolgoliubov coefficients turn out to be different. 



In pj|. Hawking radiation is obtained semiclassically by considering the quantization 
of scalar field modes in the background of a BTZ spacetime. The modes are shown to 
have the asymptotic form which resembles the standard plane wave form with a decay 
factor. The behavior at spatial infinity had to be taken care of by using reflecting boundary 
conditions since the AdS spatial infinity is timelike. The scalar field is equated to zero at the 
spatial infinity. In our case, we have the wave functionals for dust. In order to evaluate the 
Bogoliubov coefficients, we assumed a scalar product of the form given in (11011) . The integral 
is evaluated over the spatial slice of constant Killing time. The value of the wave functional 
for each shell is asymptotically decaying for large R. This makes the wave functional zero 
at spatial boundary of the AdS. Thus our solution automatically addresses the problem of 
the AdS spatial infinity being timelike. 

There are several outstanding issues that must remain for future work. One of particular 
interest is to obtain a microscopic description of the entropy of the BTZ black hole. The 
BTZ black hole is a special case of our solutions, corresponding to a constant mass function. 
Its entropy was computed using the AdS/CFT correspondence, although the physical sig- 
nificance of the degrees of freedom being counted in that approach remains ambiguous. An 
advantage of the canonical approach is that the degrees of freedom have a transparent phys- 
ical meaning. A comparison between the microscopic degrees of freedom from the canonical 
theory and those counted using the AdS / CFT would be illuminating. 
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Another issue concerns the study of a possible singularity avoidance in the simplified realm 
of 2+1 dimensions. Singularity avoidance has been claimed from loop quantum gravity [sf 



and has been shown to hold for collapsing dust shells in 3+1 dimensions SOj. The quantum 



collapse of shells in 2+1 dimensions was recently examined in [31j and the results suggest 
that at least a portion of the collapsing shell rebounds, reemerging from the horizon in a 
finite proper time, while the remaining portion of the shell is lost to Hawking radiation. 

A third issue concerns the role of naked singularities in the quantum theory. This problem 
is of particular experimental interest because semi-classical arguments suggest that naked 
singularities are quantum mechanically unstable. If those arguments hold, naked singulari- 
ties could be experimental windows into the world of quantum gravity. 

Finally, the fact that the WKB approximation becomes exact in the lattice regularization 
means that it is possibly useful to look for alternative regularization schemes. We plan to 
address some of these issues in future publications. 
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APPENDIX : Uniqueness of solutions 

In this appendix we discuss tlie uniqueness of the solutions in ( |95|) . We will consider the 
solutions in the exterior (JF > 0). The solutions in the interior are obtained in an analogous 
manner by replacing the trigonometric substitution below by the analogous hyperbolic one. 

For JF > 0, the Hamilton- Jacobi equation [the first of ([H2D] is solved by taking 



dWj cos rjj 

(118) 



dWj sin rjj 



dRj J^j 
Integrability of this solution implies that 



8r, -V'^''*"*4-1^- 
but, inserting (11181) into the second equation of (1821) . we find another equation for rj^r, R, F): 
siinij drij drij 91njF,- Aj , . 



which, using the relationship between A and the measure derived earlier, simplifies to 

^ _ ^/^eotr?.^ + ./^^^ (121) 

However, equations (I12ip and (I119P are consistent if and only if 

tanr^, = (122) 

where = aj(Tj) is a function only of tj. Taking derivatives with respect to tj and Rj 
respectively and reinserting them into either the integrability condition (I119p or into (I12ip 
now gives an equation for a(r), 

dTj fi] dRj 2 dRj ■ ^ ^ 

Now as a is a function of r but not of R, and fi, T are functions of R but not of r, the above 
equation requires a. = const, and therefore 



/i] dRj 2 dRj 
Thus we obtain the solutions (cf. fl9T|) and fl92l) ) 



(124) 



(3- f y^^'^j-^j 

Hj = — =, Wj = const. ± QjTj ± / dRj— , (125) 

1 - a^JF,. J 
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where aj and (3j — ajaj are constants. They are unique under the given conditions. 
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